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Abstract — The purpose is to calculate the heat transfer 
surface area of the shell and tube heat exchanger aiming 
recovery to waste heat with Kalina cycles. The procedure 
for calculating the area was implemented in software 
Engineering Equation Solver (EES), where the 
calculation of the LMTD was held with the existing 
functions in EES for the correction factor calculation, 
and additional correlations were implemented in the EES 
to calculate the thermodynamic properties of NH 3 -H 2 O 
mixture and to determine the transport properties of the 
phase change regions. As input data for the calculations 
were employed geometry (internal and external diameters 
of the tubes, the exchanger arrangement, geometrical 
arrangement of the tubes) and the heat balance of the 
heat exchanger (temperature, concentration ofNEJs-H20 
mixture, etc.) obtained by simulation of the Kalina cycle. 
As an example of the results that were obtained for the 
evaporator, the overall heat transfer coefficient was 
found to be 108,5 W/m 2 K and the heat transfer surface 
area of 1170 m 2 . This value of the overall heat transfer 
coefficient is within the ranges reported in the literature 
by other researchers. 

Keywords—heat exchanger, waste heat recovery, Kalina 
cycle. 

I. INTRODUCTION 

The waste heat from production processes are usually at 
low temperature. Hence, it is necessary a good heat 
exchanger design for recovery and produce power in 
thermodynamic cycles. The Kalina cycle [1-4] is 
excellent choice for waste heat recovery from production 
process. This cycle has as working fluid an ammonia- 
water (NH 3 -H 2 O) mixture that offers a better efficient on 
other cycles such as Rankine cycle. The use of NH 3 -H 2 O 
mixture is a complicating factor in determining the 
convective coefficient of heat transfer, the calculation of 
which has a certain complexity. Empirical correlations are 
used to find the convective coefficients the cold and hot 


side, and so determine the overall heat transfer 
coefficient. The determination of the heat transfer surface 
area is linked to the overall coefficient of heat transfer 
intrinsically. In this context the purpose in the design of 
shell and tube heat exchangers is to calculate the heat 
transfer surface area implemented in software 
Engineering Equation Solver (EES), aiming recovery to 
waste heat with Kalina cycles for power generation from 
waste heat of cement production process. 

II. LITERATURE REVIEW AND PROPOSED 
CYCLE 

The determination of the heat transfer coefficient for 
mixtures is complex, particularly when the process 
involves phase change. Rivera and Best [5] found 
experimentally that the heat transfer coefficient in boiling 
NH3-H2O mixture is 2 to 3 times higher than that of the 
ammonia-lithium nitrate mixture. Khir et al [6] reported 
that heat and mass flow strongly influences the heat 
transfer coefficient in the NH3-H2O mixture, but the 
concentration of the mixture does not have great 
influence. Araujo et al [7] studied the behavior of the 
coefficient of heat transfer of various fluids in the 
evaporation process, and fluid from the studies, R717 
(ammonia) presented a more efficient heat transfer. Shah 

[8] modified correlation shape to be applied in mixtures, 
and validated using the comparison with other 
correlations in literature. 

The NH3-H2O mixture has a boiling point below than 
pure water and boiling temperature variable promoting 
the reduction of losses in heat transfer, according Mirolli 

[9] , this provides an efficient use for waste heat recovery 
from exhaust gases in cement process. Kalina et al [10] 
found a large ratio between the exchange surface of the 
evaporator and the power generated by the cycle. And in 
thermodynamic cycle for waste heat recovery is known 
that the heat exchangers are devices that influence over 
the initial investment cost and the amount of power 
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generated by the system, hence the importance of a proper 
sizing of this equipment second Arrieta et al. [11]. 

The Kalina cycle proposed by Arrieta et al. [11]. for the 
recycling of waste gases from a fictitious plant cement 
production is presented in Fig. 1. The thermodynamic 
calculation results of this cycle were used as a basis for 
the development of the employed method of calculation 
the determination of the convective heat transfer 
coefficient for heat exchangers which phase shift 
occurring NH 3 -H 2 O mixture, boiling process in the 
evaporator and condensation process in the condenser. 



Fig. 1: Kaline cycle used for heat exchanger sizing 


In Fig. 1 each number corresponds to an equipment, 
being: 

- 1 turbine, 

- 2 superheater, 

- 3 evaporator, 

- 4 economizer, 

- 5 high temperature regenerator, 

- 6 low temperature regenerator, 

- 7 pump, 

- 8 separator, 

- 9 valve, 

-11 condenser. 

III. INPUT DATA IN EES 

The calculations in EES were started with the input data 
related to thermodynamic and transport properties of each 


state involved presented in Table 1 and Table 2. This data 
come from the thermodynamic calculation of the cycle 
performed by Arrieta et al. [11]. 


Table 1: Input data for evaporator 



Property 

Evaporator 


Inlet 

Outlet 


Temperaure (K) 

557. 76 

401.94 

Residual 

Pressure (kPa) 

101.32 

101.32 

Gas 


Mass flow 
(Kg/s) 

47.747 

47.747 


Fluid 

379.68 

431.45 

NH3-H2O 

Pressure (kPa) 

6375 

6375 

mixture 


Mass flow 
(Kg/s) 

0.8677 

0.8677 



Table 2: Input data fo 

r condenser 


Fluid 

Property — 

Condenser 

Inlet Outlet 


Temperaure (K) 

295.15 

303.15 

W ater 

Pressure (kPa) 

250.00 

250.00 

Mass flow 
(Kg/s) 

265.034 

265.034 


Temperaure (K) 

340.28 

300.15 

NH 3 -H 2 O 

Pressure (kPa) 

902.62 

902.62 

mixture 

Mass flow 

(Kg/s) 

0.8677 

0.8677 


IV. INTERNAL HEAT TRANSFER COEFFICIENT 

The methodology for calculating the convective transfer 
coefficient for the process of change of phase mixtures, 
both boiling as condensation, was presented by Shah [8] 
and is based on the addition of the two parameters, the 
nucleate boiling of the heat transfer mechanisms and 
effects of the forced convection. 


h mix h nb + hpc 

h NB =h l F TS (230Bo 05 ) 

Where: 
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Where in equations 1 to 4 the hmx is the convective heat 
transfer coefficient of the mixture; hm is convective heat 
transfer coefficient of the nucleate boiling; Hfc is 
convective heat transfer coefficient of the forced 
convection; hi is convective heat transfer coefficient of 
the liquid phase; Fts correction factor of the Thome and 
Shakir; Bo is boiling number; Re is Reynolds number; Pr 
is Prandtl number; ki is thermal conductivity of the liquid 
phase; Din is inside diameter; hpB is heat transfer 
coefficient of the pool boiling; q is heat flux, Tdp is dew 
point temperature of mixture; Tbp is bubble point 
temperature of mixture; pi is density of the liquid phase; 
hig is latent heat of vaporization. 

The second installment of the convective heat transfer 
coefficient of the mixture on the forced convection is: 


h = 

iipc 
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Where in equations 5 to 7 the Co is Convection number; 
Y is factor in Bell-Ghaly method; h v is convective heat 
transfer coefficient of the vapor phase; x is vapor quality; 
pi is density of the vapor phase; C P , V specific heat of vapor 
at constant pressure; H is specific enthalpy. 


V. EXTERNAL IIEM TRANSFER COEFFICIENT 

The staggered arrangement of circular tubes with fins is 
applied to the evaporator. Thulukkanam[12] proposes the 
following correlation for external heat transfer 
coefficient: 


=0.13<Re„r'‘(Pr)‘' 3 
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For staggered and smooth tubing used in the condenser, 
Bejan [13] presents a number of correlations for the 
definition of convective heat transfer coefficient for the 
outside dependent on the Reynolds number. For Reynolds 
number between 1000 and 20,000, Nusselt number is 
calculated as: 


Nud = 0.35(Re o )°' 6 (Pr)°’ 36 


Pr 

Pr 


y/4 


(9) 


Where in equations 8 and 9 the h ou t is convective heat 
transfer coefficient of the outside; s a is spacing between 
fins; h a is fin height; S a is thickness fin; Nud is Nusselt 
number; Pr w is Prandtl number in the wall temperature. 


VL CALCULATION OF THE ARIA OF THE 
HEAT TRANSFER 

The LMTD method was applied for specification of the 
thermal exchange area to each heat exchanger. The 
correction factor is calculated using proper function of the 
ESS software. 


d = U a AA T lm F 


( 10 ) 


The overall heat transfer coefficient ( U ,,) for smooth or 
finned tubes was used the formulation of Incropera et al. 
[14], which was modified to meet the particularities of the 
exchangers studied. 
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Where in equation 11 theR f /, £ is resistance to deposition of 
cold fluid; R« is thermal resistance of the wall, A c is area 
of the cold side; Rd.h is resistance to deposition of hot 
fluid; >] a is Fin efficiency; Ai, is area of the hot side. 


And AT [m given by: 
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Where in equations 10 to 14 the q is heat rate; U 0 is 
overall heat transfer coefficient; A is area of heat transfer; 
A Tim is mean logarithmic temperature difference; F is 
factor correction; Th.in is the inlet temperature of the hot 
fluid; Tc.out is the outlet temperature of the cold fluid; Th.out 
is the outlet temperature of the hot fluid; T c .i n is the intlet 
temperature of the cold fluid. 


VII. RESULTS 

The overall heat transfer coefficient is found within the 
range specified by Walas [16], for condenser type shell 
and tube of the overall coefficient of heat transfer is 
between 283.5 and 1134.9 W/m 2 K. hi the evaporator, the 
specified range is 141.75 to 340.47 W/m 2 K. Table 3 and 
Table 4 shows the values found for internal heat transfer 
coefficient, external heat transfer coefficient, overall heat 
transfer coefficient and area of the heat transfer for the 
evaporator and condenser, respectively. 


Table 3: Results for Evaporator 

Parametro Values 

Internal heat transfer coefficient (W/m 2 K) 584.8 

External heat transfer coefficient (W/m 2 K) 681.5 

Overall heat transfer coefficient (W/m 2 K) 108.4 
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Area of the heat transfer (m 2 ) 1170.0 


Table 4: Results for Condenser 

Parametro Values 

Internal heat transfer coefficient (W/m 2 K) 943.9 

External heat transfer coefficient (W/m 2 K) 3538.0 

Overall heat transfer coefficient (W/m 2 K) 660.3 

Area of the heat transfer (m 2 ) 889.2 


The parametric studies performed in EES show the 
influence on global heat transfer coefficient, heat transfer 
surface area and average convective heat transfer 
coefficient inside and outside the tubes with the variation 
of the outer diameter (0.0761 < D ou t < 0.108 m). Note 
that the outer diameter was connected such to the inner 
diameter so as to keep constant the thickness of the tube 
wall, that is the variation in the outer diameter of the inner 
diameter also varies. 

VIII. CONCLUSION 

The overall heat transfer coefficient and heat transfer area 
to the evaporator were found 108.4 m 2 and 1170.0 
W/m 2 K respectively. In the condenser, the overall heat 
transfer coefficient and area were found to 660.3 W/m 2 K 
and 889.2 m 2 respectively. For the parametric analysis, it 
is concluded that the external diameter of the tube has a 
strong influence on the overall heat transfer coefficient, hi 
the evaporator, the increased outer diameter provides an 
increase in the overall coefficient of heat transfer, since 
the reverse occurs in the condenser, mainly caused by 
constmction differences of the two heat exchangers. 

It is suggested as a future study the optimization of the 
evaporator and condenser that are used in the Kalina 
cycle, as in the case study found in Saldanha et al. [16] 
for shell and tube heat exchanger with single phase flow. 
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